Dynamic facilitation explains democratic particle motion of metabasin transitions 
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Transitions between metabasins in supercooled liquids seem to occur through rapid "democratic" 
collective particle rearrangements. Here we show that this apparent homogeneous particle motion 
is a direct consequence of dynamic facilitation. We do so by studying metabasin transitions in 
facilitated spin models and constrained lattice gases. We find that metabasin transitions occur 
through a sequence of locally facilitated events taking place over a relatively short time frame. 
When observed on small enough spatial windows these events appear sudden and homogeneous. 
Our results indicate that metabasin transitions are essentially "non-democratic" in origin and yet 
another manifestation of dynamical heterogeneity in glass formers. 



INTRODUCTION 

In this paper we study metabasin (MB) [1] transitions, 
that is, transitions between low energy/low activity con- 
figurations, in kinetically constrained models (KCMs) of 
glass formers [2-6]. We follow closely the recent work of 
Appignanesi et al. [7] (see also [8, 9]) who studied this 
problem in an atomistic model, a Lennard- Jones binary 
mixture [10], by means of molecular dynamics simula- 
tions. In Ref. [7] it was found that transitions between 
metabasins involved relatively fast and collective rear- 
rangements of a significant number of particles forming 
compact clusters in space. These apparent homogeneous 
relaxation events were termed "democratic" [7] . A natu- 
ral question is whether democratic events are distinct in 
nature to those associated with dynamic heterogeneity 
[11]. Here we show that the apparent democratic par- 
ticle motion that occurs during a metabasin transition 
is a direct result of dynamic facilitation, an intrinsically 
"non-democratic" property of KCMs which is the origin 
of dynamic heterogeneity [12, 13] in these systems. 



METABASIN TRANSITIONS IN FACILITATED 
SPIN MODELS 

The simplest microscopic models built on the idea 
of dynamic facilitation are the so-called facilitated spin 
models (FSMs), such as the Fredrickson- Andersen (FA) 
model [2] and its directional counterpart the East model 
[3]. In one dimension these models are described by a 
chain of Ising spins rti = {0; 1}, with a trivial Hamilto- 
nian H = ^ n^. Here m = 1 represents a mobile site, or 
excitation, conversely rti = represents an immobile, or 
jammed site. Glassiness is the result of local dynamical 
rules that specify the ability of a site to change state. In 
the FA model a site can only flip if either of its nearest 
neighbours are in the excited state. The transition rates 
are: 10 11, 11 -U 10, 01 11, 11 -U 01, where 
e = e _/3 and (3 = 1/T. For the East model the last two 
transitions are not allowed: a spin may only flip if its 



nearest neighbour to the left is excited, consequently ex- 
citations propagate in an eastward direction. Relaxation 
is Arrhenius in the FA and super- Arrhenius in the East 
model [2, 3, 6, 13, 14]. 

The top panels of Fig. 1 show trajectories for the FA 
and East models. Mobile sites are dark and immobile 
sites are white. In both cases the trajectory is several 
times the length of the structural relaxation time, Tq,, 
at the corresponding temperature. The trajectories il- 
lustrate the familiar features of dynamic facilitation [13]: 
at low temperatures trajectories are spatially heteroge- 
neous, excitations form continuous lines in space-time, 
and there are large inactive space-time "bubbles" . 

In analogy with Ref. [7] we study metabasin transitions 
by computing the "distance matrix" (DM), defined as: 



AHt\t") = ±Y / \k i (t>)-k i (t")\ 2 



(i) 



where ki(t) is the number of times site i has changed 
state in time t. Hence, A 2 (t',t") represents the average 
squared number of configuration changes, or "kinks", in 
the time interval (£', t"). Our definition of the DM should 
be compared to that used in Refs. [7-9, 15] for atomistic 
models, 



A V,0 = ^f>(i')-r;(0| 2 



(2) 



where r^(£) is the position of particle i at time t [16]. 
In this case K 2 (t',t") measures the system's averaged 
squared displacement between times t' and t" . 

Just as in the atomistic case, in order to focus on indi- 
vidual metabasin transitions we need to consider a small 
enough subsystem [17]. The boxes in the top panels of 
Fig. 1 indicate the sub-trajectories we use for the anal- 
ysis. The spatial extension of our subsystems in each 
case is comparable to the dynamic correlation length, 



( n iY 



e 13 [13, 17]. Due to the heterogeneous nature 



of the dynamics, mobility within a sub-region changes 
dramatically over the course of time. For long periods 
the system can remain inactive until eventually a sudden 
burst of activity occurs as an excitation line enters and 
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FIG. 1: (Top) Trajectories for the FA model (left) at inverse temperature (3 = 3, and for the East model (right) at f3 — 4. 
The length of the trajectories is several relaxation times. The boxes indicate the sub- trajectories used to calculate the distance 
matrices. The subsystem size is N — 50 for both models. (Bottom) Corresponding distance matrices A 2 (t / ,t // ) for the 
sub- trajectories in the FA model (left) the East model (right). 



eventually passes through the observation window. The 
bottom panels of Fig. 1 show the DMs for the chosen 
sub-trajectories. 

The DMs for the simple FA and East models show re- 
markable similarity to those seen in atomistic models in 
Ref. [7]. The large black squares in the DMs of Fig. 1 
correspond to times when the subsystem is trapped in a 
metabasin. It is clear from inspection of the associated 
sub-trajectories that these periods correspond to the in- 
active bubbles in space-time. These quiescent periods 
are punctuated by bursts of activity that rearrange the 
subsystem into a different metabasin, seen as new black 
square in the DM. Since our simple microscopic models 
are built purely from the idea of facilitation it is clear 
that the burst of mobility during a MB transition is not 
the result of a democratic motion of particles but rather 
a sequence of cooperative events that occur on a time 
scale far smaller than the typical MB lifetime. 

To further illustrate this idea consider the measure 
J 2 (£, #), the averaged squared kinks within the sub-region 
in a time interval 



5 2 (t,0)=A 2 (t-0/2,t + 0/2), 



(3) 



which is shown in Fig. 2 (solid lines). This quantity is 
simply A(t', t") measured along the diagonal t" = t ,J r0. 
Again, following [7] we specify that 9 is significantly 
smaller than r a but larger than the timescale associ- 
ated with any microscopic motion. For the FA model 
we choose 6 = 800 and for the East model 6 = 8 x 10 5 . 
From the figure we see that S 2 (t,6) is low during peri- 
ods of inactivity in the trajectory and there are sharp 
peaks signifying the transition between metabasins as an 
excitation line sweeps through the observed region. 

The notion of democratic particle motion [7] was in- 
ferred from the fact that the peaks in the average mo- 
bility did not result from a substantial displacement of a 
few particles, but rather an increase in the average dis- 
placement of a significant fraction (as much as 30%) of 
particles within the subsystem. Something similar occurs 
in the FSMs. In Fig. 2 we also show m(t, cV), the fraction 
of sites in the sub-region that have experienced a higher 
than average number of kinks in the time interval (j). The 
peaks in ra(t, <p) coincide with those of 5 2 (t, 0) and at the 
corresponding time a significant fraction of sites within 
the sub-region exhibit a high number of kinks. How- 
ever, this apparent democratic motion is the of result of 
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FIG. 2: Average squared kinks 5 2 (t,Q) (solid line) for the FA model (left) and East model (right). — 800 for the FA model, 
and — 8 x 10 5 for the East model. Also shown is the fraction of sites m(t, with a higher than average number of kinks 
in the time interval <j> = 200 (FA) and <j> = 2 x 10 5 (East). Dashed lines indicate the concentration of excitations within the 
sub-region, c su b(t); the curves are unsealed and have been shifted vertically to aid clarity. 



facilitation, the global relaxation event resulting from a 
sequence of locally facilitated ones. The extent to which 
the excitation line penetrates the sub-region determines 
the scale of the apparent democracy. 

As illustrated in the space-time trajectories of Fig. 1 
for both the FA and East models it is the presence of 
excitations that drives the dynamics within a given re- 
gion of the system. The features observed in the DMs 
are a direct result of the fluctuating concentration of ex- 
citations within the sub-region, c su b. To illustrate this 
point we plot c sn ^(t) along each trajectory, shown as the 
dashed line in both panels of Fig. 2 (note that the curves 
for c su b are unsealed and have been shifted vertically to 
aid visualisation). We can see that the quantities c su b(t) 
and 5 2 (t,6) behave similarly. 

For larger subsystems the sharp features seen in the 
DMs begin to disappear as the dynamics becomes more 
homogeneous in nature. Fig. 3 shows DMs for the 
East model at two larger system sizes; as the size of the 
sub-region is increased the island structure of the DMs 
is quickly lost forming a dark band along the diagonal 
t" = t' . The increasing homogeneity in the temporal 
variation of the subsystem dynamics is also reflected in 
the average squared kinks, right panel of Fig. 3. Here we 
show the average squared kinks for the East model using 
a subsystem of size N = 300, where = 8 x 10 5 once 
again. Although small in comparison to Fig. 2 the dy- 
namics within the sub-region still fluctuates significantly, 
a detail which is somewhat obscured in the DM itself. In 
addition fluctuations in the fraction of sites experiencing 
above average kinks, m(t, aV), still coincides with 5 2 (t, 0). 
The dashed line again indicates the concentration of ex- 
citations within the sub-region, the fluctuations of which 
agree precisely with ra(£, (j>). 



METABASIN TRANSITIONS IN CONSTRAINED 
LATTICE GASES 

Another class of simple KCMs are constrained lattice 
gases [4-6, 18]. Here we consider the two- vacancy as- 
sisted triangular lattice gas, or (2)-TLG [3, 19]. The 
system consists of hard-core particles that move on a 
two-dimensional lattice of triangular geometry; there are 
no static correlations between particles and at most each 
site can hold one particle at a time. Any particle on the 
lattice can only move to one of its six nearest neighbour 
sites if the following rules are satisfied: (i) the target 
site is unoccupied and (ii) both the two mutual nearest 
neighbours of the initial and target site are also empty. 
The physical interpretation of the dynamical rule is the 
steric constraint on particle motion that occurs within a 
dense fluid. The dynamics of the model is highly collec- 
tive, for a particle to be able to move first requires the 
cooperative rearrangement of many of its neighbours. 

We again choose to examine the dynamics within a 
small sub-region of a large system using the definition of 
the distance matrix given in Eq. (1). This definition is 
convenient since it accounts for all particles which enter 
or leave the sub-region during the course of a given tra- 
jectory, rather than those simply present at a given time. 
Here a kink corresponds to a particle entering or leaving 
a site. By using a sub- system within an extensive lattice 
we also avoid the potential problem of forming a back- 
bone of frozen particles that could percolate the region 
of interest [4, 5]. 

Fig. 4 shows an example distance matrix for a 10 x 10 
sub-region taken from a (2)-TLG lattice at particle den- 
sity p = 0.79. The total simulation length is 10 6 , the per- 
sistence time at the chosen density [19]. The DM exhibits 
the same features found in the FSMs and in atomistic 
models. The trajectory is characterised by extended re- 
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FIG. 3: (Left and centre) Distance matrices for the East model at two different subsystem sizes, N — 150 (left) N — 300. 
(centre). (Right) Average squared kinks 5 2 (t,6) (solid line) for the East model subsystem of size N = 300 where = 8 x 10 5 . 
Also shown is the fraction of sites m(t,<fi) with a higher than average number of kinks in the time interval <j) = 2 x 10 5 . The 
dashed line (which has been raised vertically and is unsealed) indicates the concentration of excitations within the sub-region, 
c S ub(t). The fluctuations coincide perfectly with m(t, <fi). 



gions of low activity, which appear as dark squares in the 
DM, with relatively quick transitions between them [20]. 
Fig. 5 (left) shows the average squared kinks, 5 2 (t, Q), for 
the same trajectory, where = 4 x 10 4 . Also shown is the 
fraction of sites with above average kinks in a time win- 
dow <j> = 10 4 . As for the FSMs, at the peaks in 5 2 (t,6) 
a significant fraction of sites in the sub-region are expe- 
riencing a high level of activity. From direct analysis of 
the particle trajectories we find the MB-like regions cor- 
respond to configurations in which the majority of parti- 
cles within the sub-region are blocked and remain frozen 
in position for large periods of time. A burst of motion 
results from a sequence of unlocking events which enables 
the particles to quickly rearrange themselves. After the 
effective mobility excitation passes the particles are left 
in a completely or partially frozen configuration. This 
dynamical picture is in direct agreement with previous 
analysis of the TLG models [5, 19, 21]. 

Due to the particle nature of the TLG model it is also 
possible to consider the original definition of the DM as 
given in Eq. (2), i.e. the subsystem's averaged squared 
displacement between times t f and t" . Since the problem 
of a backbone prevents simulation of a small system it is 
necessary to modify the summation such that only those 
particles present in the sub-region at time t' are consid- 
ered. Fig. 6 shows (left panel) the resulting DM for the 
same 10 x 10 subsystem used previously. The DM is sim- 
ilar to that of Fig. 4 indicating that both definitions of 
the distance matrix capture the same dynamical infor- 
mation. This is also evident when analysing the average 
squared displacement 5 2 (t,9) (again using = 4 x 10 4 ) 
as shown in the right hand panel (solid line). Also shown 
is the fraction of particles with above average displace- 
ment in the time interval <fi = 10 4 . In agreement with the 
results of atomistic simulations the peaks in the average 
squared displacement are accompanied by a large dis- 



placement of a significant fraction of the particles within 
the subsystem. 

In a recent paper [21] we identify a possible elementary 
excitation for the (2)-TLG by studying the dynamics of 
the model with the iso-configurational (IC) method [22]. 
Here one constructs an ensemble of equal length trajec- 
tories which share a common initial particle configura- 
tion but with different initial particle momenta (for the 
lattice gas Monte Carlo used to simulate the TLG mod- 
els one uses a different sequence of attempted moves for 
each trajectory). As for the Lennard- Jones liquids stud- 
ied previously [22], we found [21] heterogeneity in the 
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FIG. 4: Distance matrix A 2 (t ; ,t ;/ ) for the (2)-TLG. The data 
corresponds to a 10 x 10 subregion from a lattice at particle 
density p — 0.79 (see Fig. 5). 
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FIG. 5: (Left) Average squared kinks £ 2 (t, 6) (solid line) for the (2)-TLG model. Also shown is the fraction of sites ra(£, 0) with 
a higher than average number of kinks in the time interval <j> = 10 4 . (Right) The top figure shows the sub-region considered. 
The dynamics within this sub-system make up the distance matrix shown in Fig. 3 and in the left panel. The numbered images 
are snapshots of the particle field within the sub-region at the times indicated in the left panel. Particles which have moved 
position in the time interval (j) — 10 4 are coloured black. 



spatial distribution of the dynamic propensity, (Arf ) IC , 
the squared displacement of particle i averaged over all 
trajectories within the IC ensemble. This heterogeneity 
was found to correlate well with a structural measure, 
the connectedness of particles to extended clusters of va- 
cancies [21]. This was defined as the total sum of the 
vacancy cluster sizes to which a particle is connected to 
through its nearest neighbours [21]. Cluster connectivity 
was also found to provide a good prediction of the instan- 
taneous particle mobility within a single trajectory. The 
most mobile particles at any given time are on average 
the most connected, typically forming a ring around a 
large cluster of vacancies. To illustrate the role of clus- 
ter connectivity for the dynamics of the (2)-TLG we plot 
the average cluster connectivity of particles within the 
sub-region along the trajectory, shown as a dashed line 
in the right hand panel of Fig. 6. Once again the curve is 
unsealed and has been shifted vertically for clarity. The 
fluctuations in the cluster connectivity match well (in 
particular the highest peaks) with those of the average 
squared displacement (or average squared kinks) further 
illustrating the role of high connectivity particles as the 
relevant dynamic excitations for the (2)-TLG model. 



CONCLUSIONS 

We have shown that transitions between metabasins in 
glass formers can be understood in terms of dynamic fa- 
cilitation. The apparent "democratic" collective particle 
rearrangements are really a sequence of localised and fa- 
cilitated events. Observed on small enough lengthscales 
they just appear sudden and homogeneous. They are in 
essence another manifestation of dynamic heterogeneity. 



The distance matrix, Eq. (1), encodes the trajectories 
in an elegant way and the similarity between the DMs we 
observe here in KCMs, Figs. 1 and 4, and the ones mea- 
sured in molecular dynamics simulations [7-9] is remark- 
able. For KCMs we know that the observed structure 
in the DM is just a projection of the space-time correla- 
tions in trajectories, that is, of inactive bubbles bounded 
by excitations lines (itself a consequence of phase coex- 
istence between active and inactive dynamical phases in 
KCMs [23]). Our results here give further support to the 
view that atomistic models have trajectories with similar 
features. 
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